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ABSTRACT 

We report on Australia Telescope Compact Array observations of the ~10 5 yr old 
pulsar PSR B0906-49. In an image containing only off-pulse emission, we find a weak, 
slightly extended source coincident with the pulsar's position, which we argue is best 
interpreted as a pulsar wind nebula (PWN). A trail of emission extending behind the 
pulsar aligns with the major axis of the PWN, and implies that the pulsar is moving 
north-west with projected velocity ~60 km s" 1 , consistent with its scintillation speed. 
The consequent density we infer for the pulsar's environment is > 2 cm~ 3 , so that the 
PWN around PSR B0906-49 is confined mainly by the high density of its surroundings 
rather than by the pulsar's velocity. Other properties of the system such as the PWN's 
low luminosity and apparent steep spectrum, and the pulsar's large characteristic age, 
lead us to suggest that this nebula is substantially different from other radio PWNe, 
and may represent a transition between young pulsars with prominent radio PWNe 
and older pulsars for which no radio PWN has been detected. We recommend that 
further searches for radio PWNe should be made as here: at low frequencies and with 
the pulsed emission subtracted. 

Subject headings: ISM: general — pulsars: individual (PSR B0906-49) — radio 
continuum: ISM 
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1. Introduction 

The spin-down observed in most pulsars corresponds to a significant rate of energy loss, which 
manifests itself primarily in the form of a magnetized wind of relativistic particles (e.g. Rees & 
Gunn 1974). Under certain conditions the interaction between this wind and its surroundings 
is observable, in the form of a pulsar wind nebula (PWN). At radio frequencies PWNe fall into 



two basic classes, plerions and bow-shock nebulae: plerions (e.g. Helfand & Becker 1987) are the 
filled-center components of supernova remnants (SNRs) in which the pulsar wind is confined by 
the pressure of hot gas in the SNR interior, while bow-shock nebulae are confined by the ram 



pressure associated with their pulsar's high velocity (e.g. Frail & Kulkarni 1991; Frail et al. 1996). 
Both are characterized by significant levels of linear polarization, a centrally-peaked morphology 
and a flat spectral index (—0.3 < a < 0; S u oc u a ). 

While only a tiny fraction of a pulsar's spin-down energy goes into producing the radio 
emission from such nebulae, radio PWNe are valuable diagnostics of the properties of both the 
wind and of the pulsar itself (e.g. |Frail et al. 1996] ). Insight into the workings of such PWNe is 
limited by the fact that only six pulsars, all with ages <^ 10 5 yr, have associated radio PWNe ( [Frail 



fc Scharringhausen 1997| , hereafter FS97). Attempts to detect radio nebulae around older pulsars 



have so far been unsuccessful (e.g. |Cohen et al. 1983[) . 



In an attempt to increase the sample of radio PWNe, FS97 used the Very Large Array at 
8.4 GHz to search for PWNe around 35 pulsars of high spin-down luminosity and/or velocity, but 
found no new PWNe down to a surface brightness T& ~ 1.2 K. They concluded that only young, 
energetic pulsars produce observable radio nebulae, and that the properties of a pulsar's wind 
may change as the pulsar ages. However, a weak, compact or steep-spectrum nebula is difficult 
to detect using the approach of FS97, because the nebula is "hidden" by the emission from its 
associated pulsar. This has motivated us to attempt an alternative strategy for finding radio 
PWNe, whereby visibilities are recorded at high time resolution so that off-pulse images can be 



produced. The full results of this program are described elsewhere QStappers et al. 1998|) ; here we 



report on the discovery of an unusual PWN associated with PSR B0906-49. 

PSR B0906-49 (l = 270?3, b = -1?0) is a 107 ms pulsar ( p'Amico et al. 198S| ) whose 
characteristic age r c = 112 000 yr and spin-down luminosity E = 4.9 x 10 35 erg s^ 1 place it 
amongst the 5% youngest and most energetic of all pulsars. Hi absorption has yielded distances 
to the pulsar in the ranges 2.4-6.7 kpc QKoribalski et al. 1995| ) and 6.3-7.7 kpc (| aravanan et al. 



1996), while a distance estimate using the pulsar's dispersion measure (Taylor & Cordes 1993) 



is 6.6lo!g kpc. In future discussion we assign it a distance 7dj kpc. The pulse profile shows a 



strong main pulse and weaker interpulse, both of which are ~90% linearly polarized (Wu et al 



1993; Qiao et al. 1995] ) . Scintillation measurements imply a transverse velocity for the pulsar of 



50^ km s 1 ( |Johnston, Nicastro, &z Koribalski 1998 ). 
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2. Observations 

Observations of PSR B0906-49 were carried out with the Australia Telescope Compact Array 
(ATCA), and are described in Table |l|. In each observation, two continuum frequencies were 
observed simultaneously with 32 4-MHz channels across each; all four Stokes parameters were 
recorded. All observations were carried out in pulsar gating mode, in which visibilities were 
sampled at a time resolution of 3.3 ms (32 bins per period) and then folded at the apparent pulse 
period. These summed visibilities were then recorded every 40 seconds. 

The flux density scale of the observations was determined by observations of PKS B1934- 
638, while polarization and antenna gains were calibrated using observations every 45 min 
of PKS B0823-500. Data were edited and calibrated using the MIRIAD package, and then 
dedispersed at a dispersion measure of 179 pc cm~ 3 . The resultant smearing due to the finite 
channel width was in all cases no more than one time bin. At each frequency the field of interest 
was imaged, deconvolved using 10 4 iterations of the CLEAN algorithm, and then convolved with 
a Gaussian beam of dimensions corresponding to the diffraction limit. 



3. Results 

In Fig |l](a) we show a 1.3-GHz image of the region surrounding the pulsar, using all time bins 
from 1997 May and 1997 Aug observations. An adjacent Hn region ( preen et al. 1998| ) contributes 
considerable extended emission, which we discard by imaging using only u — v spacings greater 
than 1 kA. The pulsar is apparent as an unresolved source at position RA (J2000) 09 h 08 m 35 s , 
Dec (J2000) — 49°13'06", in agreement with its timing position ( D'Amico et al. 1988j ). Extracting 



the time dependence of the data at these coordinates from the visibilities gives the profile shown 
in Fig ^, clearly showing the main pulse and interpulse. Similar profiles are obtained at other 
frequencies: flux densities for the pulsar (including any coincident unpulsed emission) are given in 
Table |2[ An unrelated source, ATCA J090850-491301, can be seen 2'.5 to the east of the pulsar, 
with flux densities also listed in Table |2|. 

Fig |l](b) shows the same field as in Fig [l|(a), but generated by only using time bins in which 
there is no discernible pulsed emission. A weak source can now be seen at the position of the 
pulsar. While the pulsar's emission is well-fitted by a Gaussian corresponding to a point source, a 
fit to this weaker source shows it to be extended, with dimensions 15'.'0 x lO'.'O at position angle 
315°±10° (N through E). After correction for bandwidth depolarization, the source is less than 
~50% polarized. The position of its peak agrees with that of PSR B0906-49 to Of! 6 in RA and to 
0f!4 in Dec. 

The flux density of the off-pulse source as a function of frequency is detailed in Table ||. In 
observations at 1.7 GHz and 2.2 GHz the source was not detected, despite being simultaneously 
detected at each epoch at a lower frequency. For non-detections we quote upper limits to the flux 
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density, derived using the method of Green &; Scheuer (1992). This involves adding to the data 
simulated sources of the same dimensions as the off-pulse source, increasing the flux density until 
the simulated source is distinguishable above the noise. The upper limits quoted in Table [2] are 
conservative, being twice the values at which a simulated source could be discerned. 

A 1.3-GHz off-pulse image from 1997 May 12 data, made using all baselines, is shown in 
Fig [3|. As well as ATCA J090850-491301 and the off-pulse source, a faint trail of emission can be 
seen, extending ~3f 5 south-east from the pulsar's position, broadening with increasing distance 
from it before abruptly terminating. The apex of the trail coincides with the position of the 
pulsar, and the position angle defined by the bisection of the trail's two flanks is consistent with 
the direction along which the off-pulse source is extended. The flux density of the trail is ~7 mJy 
In 1.3-GHz data from 1997 Aug 17 observations, the trail is just visible amidst confusing thermal 
emission. The trail was not detected in 1.2-GHz observations (which were severely affected by 
interference), nor was it seen at 1.7 or 2.2 GHz (as for the off-pulse source). While the inclusion of 
short baselines produces a great deal of confusing structure, no other feature resembling the trail 
is seen in any part of the field in any observation. 



An Ha observation of the region (Buxton, Bessell, & Watson 1998) shows no emission which 



can be associated with any radio-emitting source, a result which can be attributed to a foreground 
dark cloud (|Hartley et al. 1986| ). 



4. Discussion 

We consider four possible explanations for the weak source coincident with PSR B0906-49: 
positional coincidence with an unrelated source, a compact SNR, emission from the pulsar itself, 
or a pulsar wind nebula associated with the pulsar. 

At A = 20 cm, extragalactic source counts (e.g. Windhorst et al. 1985) imply a probabibility 
~10 -5 of finding a background radio source of flux density >1 mJy within l" of the pulsar, 
allowing us to rule out a coincidental alignment. We estimate that positional coincidence with an 
unrelated Galactic source is also unlikely. 



The source under consideration may be a compact SNR, in the form of either a shell (Green 
1985) or a plerion ( Hclfand et al. 1989j ). However, a shell or plerion of age r c ~ 10 5 yr and 
radius ~ 0.3dj pc as required here results in an improbably high (> 10 8 cm -3 ) density for the 
surrounding medium. No other emission is seen in the region which might correspond to a more 
extended SNR ( Duncan et al. 1995 ; Green et al. 1998]) , and we presume it to have dissipated. 



Although we took care to create an image using only off-pulse bins, it is impossible to exclude 
the possibility that Fig |l](b) still contains some flux from the pulsar, in the form of either unpulsed 
emission or a further weak component to the pulse profile. Our main argument against such 
interpretations is that the off-pulse source appears extended. Although the source is too weak 
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to fit to in the u — v plane, its extension and position angle are quite different from that of the 
synthesized beam, and are visible at all reasonable contour levels. It is difficult to see how the 
presence of noise, sidelobes or other artifacts could artificially create such an appearance. The 
extension is not directed towards the phase center and so cannot be accounted for by bandwidth 
smearing. Nor is there evidence for significant time variability in the pulsar's flux, which might 
also mimic an extended source (and in any case would cause the pulsar itself to appear similarly 
extended). The source is evident in repeated observations, and persists when imaging various 
subsets of the off-pulse component of the pulse profile. The off-pulse source also lacks the high 
fractional linear polarization observed for the pulsar. We thus conclude that it is unlikely that the 
source represents emission from the pulsar itself. 

The remaining possibility is that this source represents synchrotron emission from a pulsar 
wind nebula associated with PSR B0906-49. The observed elongation of the PWN suggests that 
it is a bow-shock nebula, an interpretation supported by the trail's morphology and its alignment 
with the PWN's major axis. We interpret the trail as a synchrotron wakeQ extending behind the 
pulsar, whose inferred direction of motion is thus along PA 315°. This wake is analogous to similar 
radio structures seen trailing behind PSRs B1757-24 flFrail fc Kulkarni 1991| ) and B1853+01 ( [Frail 
et al. 1996| ) and behind the X-ray binary Circinus X-l ( [Stewart et al. 1993| ). 



By demanding that the pulsar has traversed the trail's length in its lifetime (which we equate 
to r c in all further discussion), we determine a velocity v > 60dj / sinz km s , where i is the 
inclination of the pulsar's direction of motion to the line of sight. The well-defined edge to the 
trail in the south-east suggests that the lack of emission beyond this point corresponds to a change 
in conditions rather than to synchrotron losses (which are most likely negligible over the pulsar's 
lifetime in any case). Given the agreement between the lower limit on v and the scintillation 
velocity (Johnston et al. 1998), we propose that the pulsar was born at the base of the trail and 
adopt a velocity v = 60vqq/ sinz km s _1 . The implied proper motion of 1.8vQod>7 1 mas yr _1 is too 
small to be measured directly, but could be constrained or ruled out through pulsar timing or 
VLBI. 

The pulsar wind flows freely from the pulsar until it encounters a reverse shock at radius r s 
(e.g. Cordes 19961 ), The location of this shock is unresolved in our data, and from our highest 
resolution data we estimate r s < W 17 dj cm. We assume that all the pulsar's spin-down luminosity 



goes into the pulsar wind; upper limits to its 7-ray flux ( [Thompson et al. 1994 ) make this a 



reasonable assumption to within a factor of 2. Thus, equating the energy density of the pulsar 
wind, E/Airrlc, with the ram pressure due to the pulsar's motion, pv 2 , we find an ambient 
(hydrogen) density n > 2(dyVQo/ sini) -2 cm -3 , significantly higher than for the general interstellar 
medium. A distance greater than 7 kpc is inconsistent with the pulsar's Hi absorption; thus only 
if the pulsar's velocity is significantly greater than 60 km s _1 can this discrepancy be reconciled. 



7 In further discussion the "PWN" corresponds to the compact off-pulse source, while the "trail" corresponds to 
the extended emission to the south-east. 
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If we adopt the lower distance limit determined by Koribalski et al. (1995) and adjust the pulsar 
velocity accordingly, the required density rises to n > 70/ sinz cm -3 . We note that the CO survey 
of Grabelsky et al. (1987) shows molecular material in the direction of PSR B0906-49, at velocities 
corresponding to the range determined for the pulsar from Hi absorption. Thus PWN B0906-49 
appears to be a bow-shock PWN in which the main contribution to the confining ram pressure 
appears to be the density of surrounding material rather than the pulsar's high velocity. While 



there are other instances of bow-shock PWNe around low velocity pulsars (see Cordes 1996 ) 



PSR B0906-49 is the only case in which the PWN has been detected at radio frequencies. 

The spectral index of PWN B0906-49 is poorly determined by the available data, but 
assuming a single power law over the observed range implies —3.5 <^ a <J —0.7. This is significantly 
steeper than for filled-center SNRs ( preen T99p and for the six other known radio PWN; we note 
that such a nebula would not have been detected by the 8.4 GHz search of FS97. For a = —1, the 
integrated radio luminosity of the PWN is Lr ~ lO 30 ^ erg s , corresponding to an efficiency 
€r = Lr/E ~ 2 x 10 _6 dy. These values of Lr and er are significantly lower than for any other 
radio PWN (see FS97). Of pulsars powering such nebulae, only PSR B1853+01 has a (slightly) 
lower E. 

The unusual spectrum of PWN B0906-49 may relate to a decrease in the efficiency of shock 



acceleration resulting from enhanced densities and reduced shock velocities ( Pineault et al. 1997 ). 
If a significant component of the radio emission is due to electrons accelerated at the shock front 
between the pulsar wind and its surroundings, the high density and low velocity we have proposed 
for this source might explain its steep spectrum and low luminosity. Alternatively, FS97 have 
argued that older, less energetic pulsars do not produce observable radio nebulae, and a simple 
explanation for both their and our observations is that the intrinsic spectrum of injected electrons 
steepens with age. Apart from PSR B1951+32 (r c = 107000 yr), B0906-49 is at least five times 
older than any other pulsar associated with a radio PWN. While PSR B1951+32 is only slightly 
younger than PSR B0906-49, the former system represents a complicated interaction between the 
pulsar and its associated SNR ( [Shull, Fesen, &: Saken 1989| ) , and it is difficult to separate the 
contribution of the pulsar wind from that of the supernova shock. 



5. Conclusion 

Using pulsar-gating observations at A = 20 cm, we have found an off-pulse, slightly extended 
source coincident with PSR B0906-49, which we interpet as a faint pulsar wind nebula. An 
associated trail implies a projected velocity for the pulsar of ~60 km s -1 along a position angle 
315°. The system is unusual in several ways: 

1. the nebula has a lower luminosity and steeper spectrum than any other radio PWN yet 
discovered; 
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2. PSR B0906-49 is older than any other pulsar known to power a radio PWN, while only 
PSR B 1853+01 has a lower spin-down luminosity; 

3. the PWN appears to be generated by a slow (~60 km s" 1 ) pulsar moving through a dense 
(> 2 cm' 3 ) medium. 

Thus PWN B0906-49 appears to be very different from other radio PWNe. To account for 
this, we propose that either the low velocity and high density inhibits shock acceleration, or that 
the particle spectrum injected by a pulsar into its PWN steepens with age. In the latter case, 
PWN B0906-49 may be in transition between the high luminosity, flat-spectrum PWNe seen 
around young pulsars, and the undetectable radio PWNe around older pulsars. We note that 
PWNe of low flux density and steep spectrum were heavily selected against in FS97's ungated, 
high frequency observations. We therefore recommend that future searches, particularly around 
intermediate-age pulsars, should be made as here: at low frequencies, and by imaging only the 
off-pulse emission. A collection of PWNe produced under a variety of conditions still needs to be 
accumulated before we can fully understand how the radio emission from a PWN is produced and 
how it depends on the properties of its pulsar. 

We thank the ATNF Time Assignment Committee for generously assigning an extra 
configuration, Mike Bessell and Michelle Buxton for providing Ha data on the region, and Dick 
Manchester for supplying the pulsar ephemeris. BMG acknowledges the support of an Australian 
Postgraduate Award, while BWS received support from an ANU PhD Scholarship and the ATNF 
Higher Degree Programme. The Australia Telescope is funded by the Commonwealth of Australia 
for operation as a National Facility operated by CSIRO. This research has made use of NASA's 
Astrophysics Data System Abstract Service and of the SIMBAD database, operated at CDS, 
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Date 


Array 


Maximum 


Time On 


Center Frequency (GHz) 




Configuration 


Baseline (m) 


Source (h) 


1 


2 


1997 May 12 


6.0B 


5969 


12 


1.344 


2.240 


1997 Aug 17 


1.5B 


1286 


12 


1.344 


2.240 


1997 Nov 06 


6.0C 


6000 


11 


1.216 


1.728 



Table 1: ATCA observations of PSR B0906-49. 



Frequency Flux Density (mJy) a Rms Noise 

(GHz) PSR B0906-49 PWN B0906-49 ATCA J090850-491301 (mJy beam" 1 ) a 

1.216 22 : 2 1.1 : 0.3 1.0 : 0.2 Oil 

1.344 25±2 1.5±0.3 1.1±0.2 0.06 

1.728 16±2 < 1.0 1.0±0.2 0.07 

2.240 14±2 < 1.0 1.0±0.2 0.06 



Table 2: Flux densities and sensitivity as a function of frequency. On the basis of simultaneous 
dual frequency observations, we compute an approximate spectral index for the pulsar of a ~ — 1. 
The differences in the pulsar's flux density between the three observations are most likely due to 
refractive scintillation. 



a l Jy = 1(T 26 W m~ 2 Hz" 1 
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Fig. 1. — 1.3-GHz images (1997 May and 1997 Aug data combined) of a region surrounding 
PSR B0906-49, with baselines shorter than 1 kA excluded: (a) all time bins; (b) off-pulse bins 
only. The rms noise in the images is limited by confusion to 0.06 mJy beam -1 (T& = 0.7 K). 
Contours are at -0.15 (broken), 0.15, 0.3, 0.45, 0.6, 0.75, 5, 10 and 15 mJy beam -1 . The resolution 
(as shown at the bottom right of each image) is 9'.' 4 x 7'/0, PA 0°. The "+" indicates the position 
of the phase center, RA (J2000) 09 h 08 m 50 s , Dec (J2000) -49°13 / 50". 
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Fig. 2. — Pulse profile of B0906-49 at 1.3 GHz. The solid, broken and dot-dashed lines represent 
total, linearly polarized and circularly polarized intensity respectively. The profile is 90% linearly 
and 3% circularly polarized, values close to those determined in previous measurements. 
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Fig. 3. — A 1.3-GHz image of off-pulse emission from 1997 May 12 data, using all baselines. 
The gray-scale range is to 0.5 mJy beam -1 , and the resolution (indicated at bottom right) is 
17"7 x 16"6, PA 5°. The rms noise is 0.12 mJy beam -1 . 



